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Abstract 

We report analytical results for the development of the viscous fingering 
instability in a cylindrical Hele-Shaw cell of radius a and thickness b. We 
derive a generalized version of Darcy's law in such cylindrical background, 
and find it recovers the usual Darcy's law for flow in flat, rectangular cells, 
with corrections of higher order in h/a. We focus our interest on the influence 
of cell's radius of curvature on the instability characteristics. Linear and 
slightly nonlinear flow regimes are studied through a mode-coupling analysis. 
Our analytical results reveal that linear growth rates and flnger competition 
are inhibited for increasingly larger radius of curvature. The absence of tip- 
splitting events in cylindrical cells is also discussed. 
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I. INTRODUCTION 



When a fluid is pushed by a less viscous one in a narrow space between two parallel 
plates (a device known as Hele-Shaw cell), the Saffman- Taylor instability arises |l|]. This 
hydrodynamic instability results in the complex evolution of the moving interface between 
the fluids, producing a wide range of patterns . Since the pioneering work of Saffman and 
Taylor these visually striking, viscous fingering patterns have been extensively studied, 
both theoretically and experimentally 0. 

One noteworthy distinction between experimental and theoretical work on the Saffman- 
Taylor problem is that the measurements usually reveal disturbances associated with friction 
at the sidewalls of the Hele-Shaw cell, while most theoretical studies sidestep this problem 
by assuming periodic boundary conditions. In principle, rigid-wall boundary conditions 
would be closer to an experimental realization of the problem. However, the use of pe- 
riodic boundary conditions, as opposed to rigid-wall ones in theoretical investigations of 
the Saffman- Taylor problem has been a matter of interest and debate during the last few 
decades 

In the early 1990's an attempt to measure the flow in a cell which approximates the 
periodic boundary conditions of the theory has been performed by Zhao and Maher ||5|. Their 
interesting experimental work considered gravity-driven viscous flow within a cylindrical 
Hele-Shaw cell (two coaxial cylinders separated by a small gap) of large radius of curvature. 
To drive the system gravitationally they allowed the fluids to form a stable flat interface, 
and then invert the cell to put the denser fluid on the top. The fluids used in ||^ had almost 
the same viscosity (low viscosity contrast), such that the process of finger competition leads 
to a nearly up-down symmetric interface. The main goal in Ref. was to directly compare 
viscous fingering fiow in a cylindrical cell with fiow in a conventional, fiat, rectangular cell, 
and verify if the rigid sidewalls could be viewed as having been replaced with periodic 
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boundary conditions. 

Comparison of the pattern evolution observed in the cyhndrical cell experiment with 
similar measurements in rectangular cells with sidewalls , shows that periodic boundary 
conditions yield few differences of results. Moreover, the cylindrical cell patterns also look 
very similar to computer simulation results obtained by Tryggvason and Aref [^], who used 
periodic boundary conditions in a fiat rectangular cell. 

The experimental results presented in suggest that, for large radius of curvature and 
low viscosity contrast, there are no significant changes made in the statistical properties of 
the Saffman- Taylor flow by friction at the boundaries. In other words, it seems that the 
boundaries affect the dynamics in much the same way whether they are provided by periodic 
boundary conditions or by the presence of physical sidewalls. Although the majority of 
flow features survive the elimination of sidewall friction, the authors in Ref. were not 
able to perform a meaningful test of growth rates in the linear regime due to experimental 
disturbances other than rigid-wall ones, related to the inversion of the cell. It is worth noting 
that sidewall effects are more prominent in cases with higher viscosity contrast P-[TT|. 



Despite the relevance and simplicity of the experimental work carried out in Ref. [Q, 
a theoretical analysis of the viscous fingering instability in a cylindrical Hele-Shaw cell 
is lacking in the literature. A theoretical study of flow in cylindrical cells, could yield 
new insigth into the possible causes that may contribute to the few differences of results 
detected in 0. In this work, we examine flow in cylindrical cells analytically. We consider 
the general case of arbitrary viscosity contrast and cell radius. In section |I|, we derive 
a generalized version of Darcy's law suitable to describe flow in cylindrical cells. This 
Darcy's law introduces a correction factor, which depends on the ratio between the cell's 
thickness and radius. This result enables us to express the differences of behavior between 



flow in rectangular and cylindrical cells in quantitative terms. In section III we investigate 
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the consequences of such differences in the flow dynamics by performing a mode-couphng 
analysis of the problem. We examine both linear and slightly nonlinear stages of evolution, 
and explicity show how linear growth rates and the dynamical process of finger competition 
are influenced by cylindrical geometry. The absence of finger tip-splitting in cylindrical cells 
is briefly discussed. Section |I^ presents our final remarks. 

II. DARCY'S LAW IN CYLINDRICAL ENCLOSURES 

In this section we present the physical system of interest and derive a generalized Darcy's 
law which is suited to bringing out the geometrical aspects related to viscous fluid flow in 
cylindrical passages. Consider two immiscible, incompressible, viscous fluids, flowing in a 
narrow gap of thickness b, between two long coaxial, thin right circular cylinders (cylindrical 
Hele-Shaw cell). The radius of curvature of the cylindrical cell is a (see Fig. 1). Denote 
the densities and viscosities of the lower and upper fluids, respectively as pi, rji and p2, r/2. 
The flows in fluids 1 and 2 are assumed to be irrotational, and between them there exists 
a surface tension a. The acceleration of gravity is represented by g, and points downward 
along the direction of cylinders' common axis. 

In order to derive a generalized version of Darcy's law, adjusted to describe flow in 
such confined cylindrical environment, it suffices to focus on a single fluid. The starting 
point of our calculation is a coordinate-free representation of the continuity equation for an 
incompressible fluid 

V-u = 0, (1) 

and the Navier-Stokes equation 

= -Vp + T^V^u + pg, (2) 

where u denotes the three-dimensional fluid velocity and p is the hydrodynamic pressure. 
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We consider cylindrical coordinates {r,(f>,z), where the z axis coincides with that of the 
two cylinders. Specializing to the case of flow in the z direction u^(r, = Uz{r,z)z, the 
continuity equation (|1]) leads to d{r Uz)/dz = 0. Consequently, the velocity can be written 
is a function of r only 

U2 = u{r) z, (3) 

where z denotes the unit vector along the z axis. 

Following the standard approach in Hele-Shaw problems, we restrict our attention to 
small velocity flows of viscous fluids, and neglect the inertial terms on the left-hand side 
of Eq. (0). Under such circumstances, we use the solution (H) to rewrite the Navier-Stokes 
equation as 



dp 1] d 

dz r dr 



d uir) 
r- 



dr 

Since the left-hand side is a function of z, and the right-hand side involves only r, Eq. 



(4) 



can be satisfied only if each side is equal to a constant of common value B. Imposing the 
no-slip boundary conditions at the solid cylindrical shells u{a) = u{a + b) = 0, we find the 
solution of the radial equation 



-(0 = ^ 



r^-Clogf-l -a^ 



(5) 



where C = [{a + 6)^ — a^]/ log (1 + b/a). In contrast to flow in usual flat, rectangular cells, 
observe that the velocity profile Eq. is not rigorously parabolic (or, Poiseuille-like) due 
to the presence of a logarithmic term. The profile is very close to a parabola for r/a <^ 1, 
but deviates from parabolic shape for larger values of r/a. 

Averaging the three-dimensional velocity u with respect to the transverse, radial direc- 
tion, defining v = (1/&) J^^'' u{r)dr, leads to the equation for the mean flow velocity of the 
fluid 



&^ F(b/a) 
V2^ 



dp 

dz 
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Equation (|^) is Darcy's law for flow in cylindrical Hele-Shaw cells. The function F{b/a) ^ 
1 + (1/30) {h/aY + 0[{h/aY] introduces a correction factor which measures the deviation 
from the fiat, rectangular case (a — oo). Equation (j^) recovers the usual Darcy's law for 
flow in flat cells with corrections of higher order in 6/a, introduced by 1 < F{b/a) < 2 . 

Darcy's law is the governing equation for Hele- Shaw-type flows 0. Therefore, quanti- 
tative comparisons between dynamical behavior of flow in rectangular and cylindrical cells 
must take into account the corrections introduced by Eq. (^, in addition to those caused by 
sidewall effects and cell inversion p . In contrast to the unavoidable disturbances mentioned 
in Ref. 0, which are hard to quantify accurately, Eq. allows precise determination of 
intrinsic, purely geometrical effects. 

Based on Eq. (|^) we can start understanding why the differences of results detected in 
Ref. were so small. In Ref. the authors used a cylindrical cell of radius of curvature 
a = 18 mm and thickness 6 = 1 mm, such that h/aK. 5.6 x 10^^. From Eq. (|^) this means a 
small correction of 0.01 % with respect to the flat cell case (a oo). In practical terms, we 
estimate that a ratio of roughly b/a ^ 10^^ would be enough in order for curvature effects 
to be considered negligible. 



III. LINEAR AND SLIGHTLY NONLINEAR DYNAMICS 

In this section we investigate the consequences of the changes introduced by the gener- 
alized Darcy's law in both linear and weakly nonlinear stages of the interface evolution. 
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We focus on two general questions: (i) on the linear level, what is the effect of Eq. (j^) on 
linear growth rates? (ii) concerning the onset of nonlinear effects, how do finger competition, 
up-down interfacial symmetry, and finger tip-splitting are influenced by Eq. (^? 

To study these issues we derive a second order mode-coupling differential equation for 
the interface perturbation amplitudes. We express the fluid-fluid interface as a Fourier series 



+00 



C(<^,^)= J2 Cn(^)exp(m^), (8) 

n=— 00 

where Cn(^) denotes the complex Fourier mode amplitudes and n=0, ±1, ±2, ... is the discrete 
azimuthal wave number. 

We exploit the irrotational flow condition to define the velocity potential = — V0j in 
fluids j = I and j = 2. Using the velocity potential, we evaluate Darcy's law for each 
of the fluids on the interface, subtract the resulting expressions from each other, and divide 
by the sum of the two fluids' viscosities to get the equation of motion 

Note the presence of the correction factor on the right-hand side of Eq. (^. To obtain 
we used the pressure boundary condition p2 — Pi = ctk at the interface 2; = where 
K = (l/a^)(d'^(/d(^'^)[l + (dC/dLpy]^'^^"^ is the interfacial curvature. The viscosity contrast 
A=(r72 — ??i)/(^2 + Vi)^ U = iP'g (p2 — Pi)/[12(?7i + ^72)] is a characteristic velocity. In 
Eq. we introduced dimensionless variables, scaling all lengths by the gap size 6, and 
all velocities by cr/12(?7i + 772). From now on we work, unless otherwise stated, with the 
dimensionless version of Eq. . 

Following steps similar to those performed in p^-p!^, we define Fourier expansions for 
the velocity potentials 



= ^jn{t) exp 



ft 

imp + (— 1)-^""^ — z 
a 



(10) 



which obey Laplace's equation and vanish as z ±00. We express (pj in terms of the 
perturbation amphtudes Cn by considering the kinematic boundary condition for flow in a 
cyhndrical cell. As in the fiat cell case, the kinematic condition n- vi 1^=^ = n-V2 \z=(; refers 
to the continuity of the normal velocity across the fiuid-fiuid interface. Substituting these 
relations into Eq. and Fourier transforming, yields the mode coupling equation of the 
Saffman- Taylor problem in a cylindrical Hele-Shaw cell 

Ck = xik) a + E Gik, k') Ck'Ck-k' + o{c!): (11) 

conveniently written in terms of the characteristic wave number k = n/a. The overdot 
denotes total time derivative, and 

X{k) = F{b/a) \k\[U -P] (12) 

is the dimensionless linear growth rate. The function 

G{k,k') = A\k\[l- sgn{kk')] (13) 

is the second order mode-coupling term, where the sgn function equals ±1 according to the 
sign of its argument. 

Based on mode coupling Eq. (|TI|) we now discuss some noteworthy features of both linear 
and weakly nonlinear regimes. Start with the linear growth rate (0): it is written as the 
product of the correction factor F{b/a) by X{k)rect = \k\[U — k"^], which incidentally, is the 
linear growth rate for the fiat, rectangular case |T^. This last observation indicates that, for 



fixed gap thickness b, there is a slight decrease of linear growth for increasingly larger radius 
of curvature a. The linear solution to Eq. (0) is purely exponential Ck^if) = Cfc(O) exp[A(fc)t], 
and introduces the correction factor F{b/a) into the linear, rectangular solution |T^. At the 



linear level, this correction becomes more and more important as time progresses and when 
b/a 1. 
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From Eq. (|T2D we can extract two relevant parameters: (i) the critical wave number 
kc = l/VU [defined by setting A(A;) = 0], beyond which all modes are linearly stable; 
and (ii) the fastest growing mode k* = k^/ VS which maximizes \{k), and dominates the 
initial dynamics of the interface. Note that k^ and k* show no dependence on the radius of 
curvature a. This behavior is illustrated in Fig. 2 that depicts the linear growth rate (|12|) as 
a function of k, in the flat cell limit (dashed curve) and for a cylindrical cell with a ^ b (solid 
curve). Note that for both cases the peak location and width of the band of unstable modes 
remain unchanged, independently of the value of a. Differences in behavior between flow in 
cylindrical and fiat cells are more pronounced around k*. These facts can be interpreted as 
follows: since k = n/a, if one increases the cell's radius of curvature a, the number of fingers 
n at the two-fluid interface is also increased in order to keep the fastest growing k constant. 

Now we turn our attention to the weakly nonlinear flow stage. We begin by discussing 
finger competition dynamics. It is well known that the viscosity contrast A has a crucial 
role in determining interfacial behavior for flow in flat rectangular cells [^,^^^,|12[. For low 
viscosity contrast {A ^ 0) the interface is nearly up-down symmetric, and increasingly larger 
asymmetry is observed for larger values of A {A ^ il)- Consequently, A has great influence 
on the dynamics of finger competition and pattern selection. In cylindrical cells, in addition 
to the parameter A, it is of interest to examine how finger competition dynamics is affected 
by radius of curvature a (or correspondingly, by the correction factor F{b/a)). 

In order to investigate finger competition we consider the infiuence of a fundamental 
mode, on the growth of its sub-harmonic. To do that we rewrite the net perturbation (^ 
in terms of cosine and sine modes, where the cosine ak = Ck + C-fc and sine bk = i (Cfc — C-fc) 
amplitudes are real-valued. Then, for consistent second order expressions, we replace the 
time derivative terms and bk on the right hand side of Eq. ([TT| ) by A(A;) ak and A(A;) bk, 
respectively. We consider a dominant fundamental of wave number kj = k*, and a sub- 
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harmonic of wave number kg = kf/2 and relatively weaker amplitude. Without loss of 
generality we may take a^^ > and bk^ = 0. Under these circumstances, we obtain the 
following equations of motion for sine and cosine sub-harmonic 

a-ks = + Aksttkj.] ttk,, (14) 

hs = - AksUkf] bk,. (15) 

Note that if A > and a^^ > 0, the fundamental accelerates the growth of the sub-harmonic 
cosine and inhibits the sub-harmonic sine mode. This causes increased variability among the 
lengths of fingers of less viscous fluid 1 penetrating more viscous fluid 2. This effect describes 
finger competition. As was in flat cells, interface asymmetry and finger competition are 
enhanced to a degree proportional to A. 

Although the second order term G{k,k') (see Eq. ([T3|) ) has no explicit dependence on 
a, the coupling between modes kf and kg makes the interface dynamics sensitive to the 
background cylindrical geometry. From equations (|14D and ([ISD we see that for a given 
A, the degree of competition and up-down asymmetry depend on the overall multiplicative 
term X{ks), which in turn is proportional to the correction factor F{b/a). To illustrate the 
combined influence of A and a on finger competition, we plot in Fig. 3 the effective growth 
rate Ag// = X{ks)[l + Aksttk^], taken from Eq. (0), as a function of viscosity contrast A. 
In Fig. 3 the dashed line corresponds to the flat cell limit, while the solid line expresses 
behavior for a cylindrical cell with a ~ 6. We see from Fig. 3 that the discrepancies between 
flow in cylindrical and flat cells are more noticeable for increasingly larger values of both A 
and b/a. Therefore, with respect to finger competition, the flow will be more stable with 
the increase of the cylindrical cell radius a. Geometrically speaking, we can say that the 
mean curvature H = (l/2a) of the cylindrical cell intensifies the competition among fingers 
in comparison with that of rectangular, planar flow. 
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We conclude this section by briefly discussing the possibility of occurrence of finger 
tip-splitting in cylindrical cells. Tip-splitting is related to the influence of a fundamental 
mode kf on the growth of its harmonic = 2k f |]12[. By rewriting Eq. ( pH] ) in terms 
of sine and cosine modes, and considering the coupling between kf and kh, we verify that 
the harmonic mode cannot be influenced by the fundamental. Therefore, at second order, 
there is no tendency for the fingers to split in cylindrical cells. Tip-splitting is absent for 
any value of the cell's radius of curvature a, including the flat cell limit a — > oo. This 
fact may be interpreted in geometric terms as follows: unlike finger competition behavior, 
which depends on, and varies with the cylinder's mean curvature H, finger tip-splitting is 



controlled by Gaussian curvature which is zero for a cylinder. In this sence, the absence 
of tip-sphtting in cylindrical cells was expected. 



IV. CONCLUSIONS AND PERSPECTIVES 

In this paper we investigated viscous flow in cylindrical Hele-Shaw cells analytically. The 
study of flow in such geometry requires modification of Darcy's law equation. A generalized 
version of Darcy's law was derived from first principles. It introduces pertinent corrections 
to usual Darcy's law in flat, rectangular geometry. We used a mode-coupling approach to 
examine the fluid-fluid interface evolution. Arbitrary viscosity contrast A and cell radius a 
have been considered. We deduced the following general results: on the linear level, there 
is an inhibition of growth for increasingly larger radius of curvature a; and, for sligthly 
nonlinear stages, we found that finger competition and interface asymmetry are enhanced 
for flow in cylindrical cells. In addition, we explained the absence of tip-splitting events in 
cylindrical cells. 

The study of geometry-related corrections F{b/a) for the case b/a > 1 could be inter- 
esting to study in the future. Within this limit, and following the lines of recent work by 
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Ruyer-Quil |T5[ and by Meignin et al. |T^, it would be of interest to study inertial correc- 



tions to the generalized Darcy's law (^), and examine gap size effects for the Saffman- Taylor 
instability in cylindrical cells. Experimental study in this direction could use, and take 
advantage of already existing, very good cylindrical Taylor-Couette cells. In addition, a 
thorough investigation of fully nonlinear flow stages in cylindrical cells, through extensive 
computer simulations, may reveal additional corrections and new dynamic behavior. Such 
numerical studies could provide a more meaningful confrontation between experiment and 
theory in cylindrical cells. 
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FIGURE CAPTIONS 



FIG. 1: Schematic configuration of viscous flow in a cylindrical Hele-Shaw cell. The dashed 
curve represents the unperturbed interface z = and the solid undulated curve depicts the 
perturbed interface z = ({ip,t). All other relevant quantities are defined in the text. 



FIG. 2: Variation of the dimensionless growth rate (12) as a function of /c, for f/ = 1: in 
the fiat cell limit 6/a — > (dashed curve), and in the cylindrical case 6/a — > 1 (solid curve). 
Note that for both situations kc = I and k* = 1/ ^/S. 



FIG. 3: Variation of the effective growth rate Ae// = X{ks)[l + Akgakj], from Eq. (0), as 
a function of A, for f/ = 1 and = 1: in the fiat cell limit b/a ^ (dashed line), and 
in the cylindrical case b/a —>■ 1 (solid line). For a given value of viscosity contrast A, finger 
competition increases for larger values of b/a. 
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